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Abstract 

We complete the computation of the 1-loop (a^) corrections to hard spectator 
scattering in non-leptonic B decays at leading power in A/m^ by evaluating the 
penguin amplitudes. This extends the knowledge of these next-to-next-to-leading- 
order contributions in the QCD factorization formula for B decays to a much wider 
class of final states, including all pseudoscalar-pseudoscalar, pseudoscalar-vector, 
and longitudinally polarized vector-vector final states, except final states with rj 
or rj' mesons. The new 1-loop correction is significant for the colour-suppressed 
amplitudes, but turns out to be strongly suppressed for the leading QCD penguin 
amplitude a^^. We provide numerical values of the phenomenological P/T and 
C/T amplitude ratios for the vrvr, irp and pp final states, and discuss corrections 
to several relations between electroweak penguin and tree amplitudes. 
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1 Introduction 



In the QCD factorization framework [1, 2] the matrix elements of the effective weak 
interaction operators relevant to charmless non-leptonic B decays B — > M1M2 take the 
(schematic) expression 

(MiM2|g.|fi) = F^^'^nO) Tl * /m,0m, + Tj" * fB<pB+ * ImAm, * fM2<pKh (1) 

at leading order in the 1/mfe expansion, where F'^^^^(O) denotes a. B ^ Mi form factor at 
= 0, /m decay constants, and (pM light-cone distribution amplitudes. The convolution 
kernels T/'^^ are short-distance, and can be expanded in a perturbation series in the strong 
coupling ag. Precise calculations of these kernels are required to make the framework 
predictive. This is particularly the case for direct CP asymmetries, since at leading order 
in the l/rrib expansion strong interaction phases arise only from these kernels, and only 
through loop diagrams. 

The kernels T/ are currently known at order [1, 2], which includes a 1-loop correc- 
tion to "naive factorization". The situation is different for T/^, which always involves the 
exchange of a hard-collinear gluon with virtuality m^A (A the strong interaction scale) 
with the spectator quark in the B meson. Due to the presence of this additional scale it 
factorizes further according to [3] 

= * J, (2) 

where the leading-order, 0{as), term is associated with the tree approximation to both 
the hard kernels Hi' and the hard-collinear kernel ("jet function") J. In a previous pa- 
per [4] we computed the 1-loop corrections to the hard spectator-scattering kernel Hi' for 
the (topological) "tree amplitudes" in two-body decays. Here we extend this computa- 
tion to the case of the (topological) penguin amplitudes except for certain flavour-singlet 
terms that contribute only when M2 is an 77 or 77' meson. Since J is also known at one 
loop [5, 6], our result completes the set of spectator-scattering kernels at C(q;^). 

The penguin amplitudes considered here provide the primary decay mechanism for 
b ^ s transitions. Their magnitudes and phases determine the size of the (direct) 
CP asymmetries in all charmless hadronic B decays, since the necessary interference of 
decay amplitudes carrying different weak and strong phases always involves a penguin 
amplitude. In this respect it is worth noting that the strong phases are confined to 
T/ at order a^. At O^a"^) a new source of strong phases appears in the spectator- 
scattering kernels T}' via the one-loop correction to Hi'. In [4] we found for the case 
of the topological tree amplitudes that this contribution is comparable in size to its 
0{as) counterpart in T/. Consequently it can change qualitatively the picture of CP 
asymmetries, and may be important both in accounting for experimental data and in 
predictions for yet unobserved asymmetries, as well as for disentangling possible new 
physics contributions from the standard model background. In the case of penguin 
amplitudes similarly large contributions could occur. ^ Spectator-scattering corrections 

^We note that numerically certain l/mh-supprcsscd, but "chirally enhanced" penguin amplitudes [1] 
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to one of the four penguin amplitudes have already been calculated in [7]. We discuss 
the difference between that calculation and ours in Section 3. 

The organization of the paper is as follows. In Section 2 we define the various flavour 
penguin amplitudes and discuss the diagram topologies that have to be calculated. We 
also set up the matching equations for the operators in the weak Hamiltonian to soft- 
coUinear effective theory (SCET) that define the full set of hard-scattering kernels Hj^. 
These kernels are in turn expressed in terms of a complete set of "primitive" 1-loop 
hard-spectator-scattering kernels. Some details of their computation and the results for 
the primitive kernels are given in Section 3. In Sections 4 and 5 we obtain the numerical 
values of the corrections to the penguin amplitudes aj(MiM2), aj^Ew(^i^2) {i = 3,4), 
and provide updated results for some phenomenologically important penguin-to-tree and 
electroweak penguin-to-tree amplitude ratios. We conclude in Section 6. 



2 Structure of the penguin kernels 
2.1 Flavour amplitudes 

Our goal is to evaluate matrix elements of the the effective weak Hamiltonian for 6 Z) 
transitions given by (see [2], where also numerical values of the Wilson coefficients Ci 
are given) 

^eff = % E KoKb I CiQ? + C2QI + E c,qA + h.c, 

V 2 p=u,c \ i=3... 10,77,89 / 

Ql = {Paba)v-~AiDbUb)v-A, Ql = iPbba)v~A{DaUb)v-A, 

Q3,5 = {Daba)v-A^{qb(lb)v^A, 
1 

Qifi = {Dbba)v-A y^(gagb)y=FA, 

3 

Q7,9 = {Daha)v-A^ 7^eq{qbqb)v±A, 

3 

Q8,10 = {Di,ha)v-A'^-eq{qaqb)v±A-, 

q ^ 

Qii = ~^Da^u{l + l^)F^''h, Qsg = -^Da,,{l + i,)G^''h, (3) 

with a, denoting color, D = d or s, and {qiq2)v^A = ^i7/i(l T I5)q2- denotes the 
electric charge of quark q in units of the positron charge e, and the sum over quarks 



are important, which are currently known to 0{as). The corresponding contributions are not the 
subject of the present paper. 
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extends over q = u,d, s, c, b. The definition of tlie dipole operators corresponds to tlie 
conventions Z-D^ = id^ + ggG^T^ + eeqA^ and a^^, = |[7/i,7!y]. Tlie effective Hamiltonian 
is understood to be renormalized in the NDR scheme as defined in [8] . 

To organize the flavour quantum numbers of the factorized matrix elements (1) we 
match TYeflf onto a transition operator such that its matrix element is given by [9] 

{M[M^\n.,\B) = Y: {M'M\n\B)- (4) 

p=u,c 

There are six different flavour structures required, resulting in 
U = 5p„ai(MiM2)A([g,M][MD]) + 5p„«2(MiM2) 

+ a§(MiM2) ^A([g,D][gg]) + af(MiM2) J] A([g,g][gZ}]) 

9 9 
9 ^ 

+ <Ew(MiM2) E|e,A([g,g][gD]), (5) 

9 ^ 

where the sums now extend only over q = u,d,s, and denotes the spectator anti-quark 
in the B meson. The coefficients af (M1M2) contain all dynamical information, while the 
arguments of A encode the flavour composition of the final state and hence determine 
the final state to which a given term can contribute. The at parameters introduced in 
[9] are in close correspondence with the widely used "graphical" or " topological" ampli- 
tudes [10]: ai (02) with the colour- allowed (colour-suppressed) tree amplitude; with 
the QCD penguin amplitude; af with the QCD flavour-singlet penguin amplitude; and 
•^SEW i.<^4Ew) with the colour- allowed (colour-suppressed) electroweak penguin ampli- 
tude. We define 

{M[M^\a^{M,M2) A{[. ..][.. .])\BJ ^ «f(M{M^) A^^m^ (6) 

whenever the quark flavours of the flrst and second square bracket match those of M[ 
and M2, respectively, and 

(M(M^|af(MiM2) A{[. ..][.. .])\B,^) ^ a^{M!,M[) Am'^m[ (7) 

whenever the quark flavours of the flrst and second square bracket match those of M'^ 
and M[. The quantity Am^M2 is given by 

(Mi,M2 = P) 
{M, = V,M^ = P) 
(Ml = P, M2 = r) 
(Mi,M2 = V) 

Here /+ and Aq denote pseudoscalar (P) and vector {V) meson form factors in the 
standard convention, and /mj are the (longitudinal) decay constants. Here and in the 
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remainder of the paper we consider only the longitudinal polarization state of the vector 
meson. This is sufficient for B PV decays, but not for B VV, where two further 
transverse amplitudes are required for a complete description. The transverse amplitudes 
are, however, 1/mfe-suppressed except for a certain electromagnetic amplitude [11]. (In 
(8) mvi{n+ ■ ^*Mi) is ^{''^b)-, hence all four expressions are of the same order in the 
heavy-quark expansion.) 

In order to exemplify the notation consider the decay Bd tt^P^, for which qs = d 
and D = d. The spectator quark can go to either one of the two mesons, so Mi can be 
7r° or p°. Hence, e.g. 

{^'p'\al{M,M2)Y^A{\dq][qd\)\Bd) = ^ laliTT^p') A^o^o + al{p%') A^oJ . (9) 



On the other hand. 



{7r'pV,{M,M,)Y.Ai[dd][qq])m = 0, (10) 



since q = u,d contribute equally but with opposite sign for the mesons vr^ and p^. We 
will assume isospin symmetry for the hadronic parameters of the mesons. It is then 
conventional to express the form factors and decays constants through one representa- 
tive member of the isospin multiplet, for instance of the charged pion in the case of 
pions. With this convention (6) must be supplied with isospin Clebsch-Gordan coeffi- 
cients 1, ±1/a/2 etc. from the flavour composition of the mesons. This is the origin of 
the factor 1/2 in (9). 



2.2 Diagram topologies 

Let us denote, in Figure 1 and below, the initial and flnal state mesons by their flavour 
content. For instance, in Figure 1(a) the B meson has the flavour quantum numbers 
[qsb], the flnal-state meson M2 represented by the up-going quark lines has flavour content 
[uD], and the remaining meson Mi [qsu], where the u line is drawn to the right. (The 
spectator line is not shown in the flgure.) 

The flavour amplitudes af in (5) receive perturbative contributions due to several 
ways of contracting the quark lines in the operators Qi with the valence quarks of the 
initial and flnal states, and possibly with each other. Considering flrst the case where 
M2 does not have a flavour- singlet component, i.e. disregarding for the moment, 
the contractions are shown in Figure 1. Here arbitrarily many gluons, that are not 
shown, may connect the quark lines, or originate from a quark line and connect with 
the spectator quark. We also consider photon exchange to lowest order, as indicated in 
Figure 1 (f) through (m). The last four contractions contribute only if a vector meson 
is involved, due to the intermediate photon. 

Focusing for the moment on the operator Qi, Figure 1(a), which we call the "right 
insertion", gives the contribution to the colour-allowed tree amplitude ai as can be 
seen by comparing the flavour labels of Figure 1(a) with (5). Contraction (b), the 
"wrong insertion", gives the contribution to the colour-suppressed tree amplitude a2, 
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Figure 1: The various contractions of the external hues with the weak effective vertex, 
excluding arbitrarily many additional gluons. The connected fermion lines indicate the 
contraction of spinor indices. Terminology: (a) "right insertion" , (b) "wrong insertion" , 
(c) [(f)] "connected [photon] penguin", (d) [(g)] "disconnected [photon] penguin", (e) 
"magnetic insertion" . The photon can attach with its 'loose' end anywhere except on the 
line it originates from. There are several more topologies (h) through (m) contributing 
to electroweak penguin amplitudes (see text), and some topologies (not shown) only 
contributing when at least one final-state meson is an SU(3) flavour- singlet. For a more 
detailed explanation, see the text. 
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while the "connected penguin" (c) contraction contributes to the (topological) penguin 
amplitude a^. For the operator Q2, the roles of the "right" and "wrong" insertions 
are interchanged. The QCD penguin operators Q3 ■ ■ - Qq contribute, through all of the 
contractions (b) (with relabeling u ^ D, u ^ q, and D ^ q) through (d), to the penguin 
amplitude 0^4 but not to ai^2- This is because of the sum over quarks already present in 
the effective weak interaction operator. The magnetic penguin operator Q^g contributes 
to a\ through contraction (e). In case of the electroweak penguin operators Q-j . . . Qio 
the right insertions (a) contribute to afEw ^"^^ wrong insertions to 04 ew However, 
insertions of Qj . . . Qiq into (c) and (d) do not contribute to the electroweak penguin 
amplitudes, but to electroweak corrections to the QCD penguin amplitudes, which we 
neglect. Thus, the second line of the figure corresponds to contributions to 04, while the 
first line can contribute to any of the af coefficients depending on the operator that is 
inserted. 

The flavour- singlet QCD penguin amplitude af is more complicated [12]. There exists 
a contribution from the right insertion (a) which is similar to the ones discussed above, 
and which will be given below. In addition, however, there are penguin contractions, in 
which meson M2 is made up only from gluons. The main reason for not calculating these 
contractions is that in case of the fiavour-singlet penguin amplitude spectator scattering 
does not factorize in the form (1), but requires the introduction of a generalized non-local 
form factor [12].^ Since the value of this form factor is unknown, the calculation of is 
rather uncertain, and it is not useful to calculate loop corrections. 

We also consider diagrams with photon exchange, but a complete calculation of QED 
corrections is beyond the scope of this paper. In fact, QED corrections to naive factor- 
ization have not even been considered for the simpler non-spectator contributions to the 
factorization formula up to now, except for an estimate of the soft photon contribu- 
tion [15]. In general, electromagnetic corrections lead to isospin violation, incorporation 
of which requires an extension of the parameterization (5). Here, similar to [2, 9, 16] 
for the non-spectator contributions, we restrict ourselves to photon exchange diagrams 
that directly give the charge and flavour structure of af ew ^^^1 a^Ew These "direct 
electroweak penguin" contributions are shown in the third and fourth line of Figure 1. 
When a photon is exchanged (besides, possibly, gluons) between the quark loop and 
the qq quark line, a contribution from Qi . . . to the electroweak penguin amplitude 
q;4Ew(-^i-^2) arises for contractions (f) and (g), and to q;|ew(^i^2) through contrac- 
tions (j) and (k), if M2 is a vector meson. Finally, the magnetic dipole operators Qr^, Qsg 
contribute to a^Ew (^'i) ^'^d, if M2 = V , to ew contractions (1) and (m). To make 
the distinction between "direct electroweak penguin" contributions and the remaining 
ones clearer we show in the first line of Figure 2 two diagrams that are included in our 
calculation of a^Ew ^^e other hand, the diagrams in the second line with the role of 

^In soft-collinear effective theory this term is related to the SCETj operator with field content 
[fW^ci^t)][W^l2*-^c2^^c2]- This operator is relevant despite the colour-octet structure of the collinear-2 
field product due to the non-decoupling of soft gluons at the level of power-suppressed interactions [13] 
in soft-collinear effective theory (SCET). The SCET rederivation of the QCD factorization formula for 
flavour-singlet mesons given in [14] is not correct, because it neglects this contribution. 
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Figure 2: Examples of closely related diagrams that are interpreted as "direct electroweak 
penguin" contributions (first line), and as isospin- violating contributions (second line). 

gluon and photon exchanged are not included, because they constitute isospin-violating 
contributions to the QCD penguin amplitude a^. 

In summary, we see that there are a number of different short-distance coefficients, 
which in general receive contributions from several weak-interaction operators Qi. The 
coefficients depend on (i) the Dirac structure of the operator Qi, (ii) its colour structure, 
(iii) the type of contraction, (iv) for the penguin and magnetic contractions, whether a 
photon is exchanged or not, and (v) on the mass of the quark in the fermion loop. 

2.3 Matching onto SCETj 

The hard-scattering kernels T/'^^ can be determined by matching QCD— *-SCETi (result- 
ing in T/, Hj^), and in the case of T^^ subsequently SCETi^SCETn (resulting in J), 
whereby the kernels are identified with Wilson coefficients multiplying non-local opera- 
tors in the effective theories [3, 4, 13, 17, 18, 19]. The two steps are equivalent to ex- 
tracting, respectively, the hard and hard-collinear momentum regions from quark decay 
amplitudes according to the strategy of expanding Feynman diagrams by regions [20]. 
Single and double logarithms of ratios of scales are summed through renormalization 
group equations in SCETj, and for the light-cone distribution amplitudes. Already in 
SCETi modes corresponding to the two different light-like directions decouple at leading 
power in A/rrib, except for the contribution to the flavour- singlet QCD penguin ampli- 
tude discussed in the previous subsection. Because of this, in SCETj there are only 
right insertions, which moreover factor into matrix elements of two currents. 

Following the notations and conventions of [13, 18], meson Mi, which picks up the 
spectator anti-quark from the B meson, moves in the direction of the light-like vector 
n_. The coUinear quark field for this direction is denoted by ^ with = 0, the 
corresponding coUinear gluon field is A^i- The second meson M2 moves in the opposite 
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direction ?7,+ , and the coUinear fields for this direction are Xi satisfying r/+x = 0, and 
Ac2- The heavy quark field is labeled by the time-like vector v = (n_ + n+)/2 with 
v"^ = 1, and i^h^ = h^. 

In [4] we argued that, ignoring colour and flavour, to leading power in l/m^, only two 
operators in SCETi are needed to match the operators Q12, as long as only current- 
current diagrams (contractions (a), (b) in Figure 1) are considered. The identification 
of these operators was based solely on power-counting arguments in SCETi [18], except 
for the fact that only left-chiral n+-collinear quark fields x were considered due to the 
{V — A) X {V — A) structure of Qi 2- Consequently, in addition to the extended flavour 
structures already discussed, the only novelty in the presence of the full set of operators 
Qi and diagram topologies is the appearance of a second leading operator in the coUinear- 
2 sector, of opposite chirality, given by {xWc2){tn_)rjl^{l + -f5){W}2X){0)- Mixed-chirality 
bilinears built from the x fields necessarily involve either additional transverse derivatives 
or a Dirac matrix carrying an uncontracted transverse Lorentz index. In the former 
case they are power suppressed, while in the latter case the resultant operator cannot 
contribute to pseudoscalar or longitudinally polarized vector mesons. Once again, the 
flavour- singlet penguin amplitude is special and requires operators with n+-collinear 
gluons. Thus, excepting af, the four SCETj operators relevant to our calculation read 



nib 



o'Ut, s) = — 



^ Wc2)(tn_)^(lT75)(H^k) 



rrib 



X 



(11) 



where 0\ = and 0^1 = O^^ in [4]. The operators 0\ ^ include the short-distance 
^f+^\^) [6] such that their matrix elements are proportional to the 



coefficients C^f^^\ 



form factor /f ^^0) (^^^'(0) for vector mesons) in QCD (not SCETi). As usual, in (11) 
fields without position argument are at x = 0, and the field products within the large 
brackets are colour-singlets. When including electromagnetic interactions the "direct 
electroweak penguins" contribute to the coefficient functions of O^]^, but there are also 
SCETi operators of the above form with the transverse collinear gluon field W^iiIp±ciWci 
replaced by a photon. This provides the sort of electromagnetic isospin-breaking, which 
as explained above we do not consider here. We can now take care of flavour by adding 
labels to the operators, 

^ C>Y}^{[q,qM^] [qM^qM^]) , (12) 

where the labels gMi,2 and qM2 give the flavours of the flelds ^, x, and respectively, and 
the redundant spectator label has been added to match the notation (5). Then up the 
n_|_-collinear gluon operators mentioned above, whose matrix elements contribute only 
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to flavour- singlet M2 (via a modification of a^), at leading power in l/mt, the complete 
weak Hamiltonian (3) can be accounted for in SCETi by 

^eff = % E y;DVr>bUu hi * Ol{[qsu][uD]) + Hi' * 0'l{[qsu][uD]) 

V p=u,c \ L 

+Tl * Ol{[qsD][uu]) + Hi' * Oii([g,D]H)| 

+ E (r3\'*Eol(M[gg]) + i^^l'^*Eol^(fc^][?"?]) 

fc=L,_R { q q 

+Tli,EW * E le.OlilqMm) + ^Sw * E le,Ol}{[qsD][qq])^^ 

+ E (rif*Eol(fcg]M) + <;^*Eol^(fc?][?~^]) (13) 

k=L,R { q q 

+tZew * E ^e,0^(feg][gD]) + if^i'V * E ^e,0^^(feg][gD])|) 

where we employed the notation 

Ti*Ol = J difl,{t)Ol{t), Hi *0'^ = J didsHllii s)0'^{t, s) (14) 

with s = n+p's = niBS, t = n^qt = rriBt, and p' (g) the momentum of Mi (M2). As in 
(5), the sums over q extend only over the light quarks u, d, s, eventually implying that we 
neglect the "intrinsic charm" content of the mesons. Of the various matching coefficients 
in (13), the 7;\(m) = J di e'""^ fl,^{i) are all known to the 1-loop order (a,) [1, 2, 9]. The 
1-loop (a^) corrections to 

Hll^{u, v) = j dids e*("*+(i-'^)^) Hl^{i, s) (15) 

have been computed in [4]. In this paper we will compute the remaining coefficients 
H'I^^ , -f^sfcEW -^IfcEW parts of -^3^'^. We do not perform an expansion in rric/mi, in 
the matching calculation. The hard matching coefficients are therefore functions of the 
ratio Sc = ml/ml, whenever diagrams with internal charm-quark loops contribute. 

The individual terms in (13) are in close correspondence with the (M1M2) ampli- 
tude parameters. The precise connection follows by evaluating the matrix element of 
(13). Because the SCET Lagrangian contains no leading-power interactions between the 
collinear-2 and coUinear-l fields, the matrix elements of 0\{t, s), Ojlit, s) fall apart into 
two factors each. For a pseudoscalar M2 = P, 

(P|(xFF,2)(tn_)^(l±75)(l^ek)|0) = T-^^ I'due^'^'Mu), (16) 
z z JO 
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while for a vector M2 = V with polarization vector we have 



such that only the longitudinal polarization state contributes. (The apparent suppression 
mv/mB is cancelled by the polarization vector.) Here 4>p{u) ((/)y(M)) denotes the leading- 
twist light-cone distribution amplitude of a pseudoscalar (longitudinally polarized vector) 
meson. Using [6] 

{P\m,,)^[Wyip^,W,^]{sn+){l + j,)K\B) = -m,mB f' dre'^'Epir), (18) 
2 Jo 

{V\i^W,i)^[Wyi^^,W,,]{sn+){l + ^,)K\B) = m,mv{n+ ■ e*) /'dre^" """^ 
I Jo 



and defining 



2m 



Sp(r) Ml = P 
SM,(r) = <{ ms ^ , . , . (19) 



2m 



r Mi = V 



V 



we obtain 



Jo 

~ 2irBA/i(o) dudvH]l'P{u,v) Sa^^(1 - v)(j)M,{u) (20) 

with F^^'^i(O) = /+*^H0) when Mi = P, and P^^^i(O) = A^*^^(0) for Mi = V^. The 
index i applies to i = 1, 2, 3, 4, SEW, 4EW, and k = L,R except for i = 1,2, where k 
is empty (as is p). The two coefficients ai(MiM2), a2(MiM2) correspond to the "tree" 
ffavour-amplitudes in (5), while the four "penguin" amplitudes are given by 

ai(MiM2) = aL(^i^2) T a^sniMiM^), 

<EW(^1^2) = «L,EW(^1^2) Tal/J,EW(^1^2), 

<EW(^1^2) = <,EW(^1^2) T a^R,Ew(^1^2). (21) 

Here the upper (lower) signs correspond to the case M2 = P(y).^ The generalized form 
factors S(r) factorize into light-cone distribution amplitudes after matching them to 



•^In the older Ui notation the coefficients 4 of the left-handed operators are related to 

a3,4,9,i0j the right-handed 3ew ^'^ '^^J- remaining a^j^, '^arew vanish. The coefficients og, ag 
in the a; notation correspond to power-suppressed penguin amplitudes. Because they have the same 
flavour structure, they are included in the definition of and a4^EW: respectively, in [9]. Similarly, 
the contributions from gluon operators omitted above arc included in a^. This will be assumed in the 
following. See also (76) below. 
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Table 1: Definition of tlie primitive kernels. Where two kernels are given, the first 
corresponds to the colour structure of Qi, the second to that of Q2- The "star" means 
that there are no corresponding kernels at leading power. Kernels in parentheses vanish 
at the 1-loop order. 

SCETn [18]. The resuh is [6] 

S(r) = - — / — fB4>B+{^) dw fM,4>MAw) J{t]w,uj) (22) 
Arrih Jo uj Jo 

with J = J||, the hard-coUinear kernel, known to order [5, 6], and J'b the static B 
meson decay constant as defined in [6]. 

2.4 Primitive kernels 

As explained at the end of Section 2.1, the short-distance coefficients encountered in 
matching an operator Qi in (3) can only depend on the Dirac structure, the colour struc- 
ture, the type of contraction, the mass of the quark in penguin loops as far as present, 
and in the case of penguin contractions and of Qsg, whether a photon is exchanged with 
the quark line on the right or not. Stripping the operators of all light-flavour and colour 
labels and an overall normalization factor, we encounter the fuU-QCD Dirac structures 

Pll = [qYil-l5)q][qi,{l-lM 
Prl = [qYil + l5)q][qi,il-l5)b], 

^mag = qcTf^uil + Ts)^, (23) 
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where q denotes some light quark field (no summation implied) and g = e ot Qg. Their 
insertions into the contractions in Figure 1 define primitive kernels according to Table 1. 

Of all these kernels, ri and r2 have been given in [4].^ The fact that these are sufficient 
to account for both "right" and "wrong" insertions (a) and (b) can be traced to Fierz 
symmetry. This is only valid because of the Fierz properties of the scheme used to define 
the renormalized effective weak Hamiltonian. In the case of penguin contractions, we will 
find, for instance, S3 7^ S2, although the differences between such naively Fierz related 
kernels always assume a very simple form. In the present paper we also calculate the 
insertion (a) of the Pjn structure. Insertion (b) results in a l/m;, correction (related to 
the infamous power-suppressed but "chirally enhanced" scalar penguin amplitudes), so 
only two primitive kernels rs 4 are needed. On the other hand, the fact that the penguin 
contractions (c) and (d) contribute for both V — AxV — A and V — AxV + A operators 
and for both colour structures requires the introduction of eight kernels Sj, s-. 



,4. 



The remaining kernels are related to electroweak penguin terms and the insertions of the 
magnetic dipole operators. 

We can now express the hard-scattering kernels Hjl in (13) in terms of these primitive 
building blocks as (n/ = n/ — 2 = 3 the number of flavours treated as massless) 



7:7-II,p 



2 



ri 



An 



+ 



Air 



CMsp) + C^s^isp) + C3(s3(0) + 53(1)) + ^4(54(0) + 34(1 
+ C54(1)+C64(1) 

+ (si(l) + si(se) + %si(0))C3 + (s;(l) + s[isc) + nfs[iO))C5 
+ (s2(l) + S2(s,) + njS2(0))C4 + (4(1) + 4(5,) + nfs'^{0))CQ 

+ miCsg 



ttU,P 
^4R 



(24) 
(25) 



o-II,P _ - 
J^3L - — 



2 



1 tts 
u Air 



C, + /r2C3 
An 



(26) 



"'The calculation of r\^2 has been repeated in [21]. There is a difference in the results for ri, which 
is simply the colour factor [Cf — Ca/'^) times the tree- level kernel. Such a difference is likely to origi- 
nate from an inconsistent treatment of ultraviolet or infrared singularities. While [4] uses dimensional 
regularization for infrared singularities, and therefore has to deal with evanescent operators, off-shell 
IR regularization is employed in [21], the price for which is the calculation of non-trivial SCETi matrix 
elements. If no error is made, the result for the hard-scattering kernels should be the same in both 
methods. In addition, in [21] a projection on the Dirac structure of the SCETi operator is taken. It is 
not obvious that such projections commute with MS rcnormalization, and in general they do not. Note 
added: After the submission of the present paper for publication, a corrected version of [21] appeared, 
which now agrees with the results of [4]. 
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H. 



U,p 
3R 



4L,EW 



H, 



4/?,EW 



H. 



ii,p 



3R,EW 



2 



1 

n An 



(27) 



2 frl 

- + 



N, 



Cg + -7—f 2C1Q 



^ 2 Clem 



e«CiS3,Ew(Sp) + e„C2S4,Ew(Sp) + edC3(^S3,Ew(0) + S3,EW 

1)) + edCsSg Ew(l) + edC6S4,Ew(l! 

C3 + (edS2,Ew(l) + e„S2,Ew(Sc 



3 47r 

+ (edSi,Ew(l) + e«si,Ew(scJ Jt-^s 

l,Ew('^c))C*5 + (crf'52,Ew(-'-) + 6«^2,EW 
(cfT 



+( (1) + 



C4 



0, 



(28) 
(29) 



2 

, 2 Ocm 





rl 


OCs 1 


{ 







Ait 



3 47r 

+ (edS8,Ew(l) + enS8,Ew(Sc 

"YCff 



/(1))C3 



e«S5,Ew('5p)Ci + erf(s5^Ew(0) + Ss^EW 

C4 + (crf'58,Ew(-'-) + ^w^8,EW 



(30) 



2 



^^3 



a. 



An 



2 a. 



+ - 



3 47r 



e«S5,Ew(Sp)Ci + erf(s5,Ew(0) + S5,ew(1))C'3 



+ (edS8,Ew(l) + euSs,Ew{Sc))C4,+ {edS8Ew(l) + e«S8Ew('5c))C'( 



+ "^3,EwC*77 + "^4,EwC'8g | • 



(31) 



Here the argument of the primitive kernels denotes the quark mass ratio Sg = m^/m^ 
with q the quark flavour in the fermion loop of the penguin contraction. We also used 
2e(i + e^j = to simplify some of the electroweak penguin kernels, and introduced the 
standard "effective" Wilson coefficients, C^^^ = Cj^ + eaiC^ + N^C^), Cfg = Csg + C5 to 
combine some ultraviolet contributions from fermion loops with the magnetic penguin 
coefficients. This convention modifies the primed primitive kernels such that Sg ^y^, Sg 

vanish at one loop. For -ffs^^EW ^srew^ terms on the last three lines of each 
that stem from the insertions (j) to (1) in Figure 1 are identical, so they cancel out for 
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M2 = P as they should. To simphfy the notation for H^l^-^y^ and H^£■^^J^^, we already 
dropped the primitive kernels that vanish at one loop. 

3 Computation and results 

In this section we provide some technical remarks on the calculation and then sum- 
marize the results for the primitive kernels. We calculate the quark-gluon amplitude 
b{p) — > q{qi)q{q2)q{Pi)g{p'2) with insertions of Qi in QCD at 1-loop and subtract the cor- 
responding SCETi matrix elements to obtain the hard-scattering coefficient. Insertions 
of the magnetic dipole operators require only a tree-level calculation. At leading order 
in the l/rrib expansion it is almost always possible to approximate the external momenta 
by their leading components p = mtv, qi = unibn+/2, q2 = umbn+/2, p[ = wmfen_/2, 
P2 = um;,n_/2, which makes the calculation particularly simple in dimensional regu- 
larization. Exceptions to this are provided by diagrams with intermediate lines whose 
virtuality vanishes in this approximation, such as (a6), (e4) and all the diagrams in the 
third line of Figure 3. 

3.1 Technical details 

Vertex contractions (Figure 1 (a),(b)) 

The calculation of the new kernels r^^^ related to insertions of {V — A)x (V + A) operators 
proceeds analogously to the one of ri^2 described in detail in [4]. Only the right insertions 
need to be computed, since the wrong insertions match to power-suppressed SCETi 
operators. Substituting 1— 75 — I+75 in Eq. (17) of [4] to deal with the new 0}j operators 
leads to a SCET operator basis, in which O2-4 are no longer evanescent (vanishing in 
four dimensions). We can remedy this by choosing a new basis, in which all 7^ stand 
to the right of all other transverse Dirac matrices. The treatment of IR singularities 
can then be done as in [4]. In the calculation of the matrix element of the evanescent 
operator O2, one now encounters in addition to the contribution shown in Figure 4 of [4] 
a contribution related to the matrix element of the ^A±cihv part of the operator. This 
contribution can be deduced from the ultraviolet pole of B in Eq. (39) of [6]. 

QCD penguin contractions (Figure 1 (c) — (e)) 

The penguin contraction diagrams are shown in Figure 3, of which the first two lines are 
relevant to the QCD penguin amplitude 04. The calculation of these diagrams is much 
simpler than the calculation of the vertex contractions. In particular, they are infrared 
finite, and no evanescent operators in SCETi need to be considered. As in the case of 
the vertex contractions, one must take into account that the QCD factorization formula 
uses full-QCD form factors by convention, hence the second line in the definition of ^ 
in (11). As a consequence, the "form-factor subtraction" 

2^i(i)^(.i)(0) ^ (32) 
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Figure 3: Summary of penguin-contraction diagrams. The 'loose' gluon or photon can 
attach to the crosses. Corresponding to diagram classes (a)-(f) there exist the diagrams 
with the horizontal gluon line replaced by a photon line. (For these (a7), (e5) and (e6) 
are absent.) 



(cf. Eq. (26) of [4]) has to be added to the diagrams in Figure 3 to obtain the final 
result for Hjl, where T^l^^ is the penguin contraction contribution to the 1-loop kernel 
Tj^. The form factor subtraction affects the structure of the result. For instance, we find 
that diagrams (e) and (f) in Figure 3 cancel exactly against the Cgg term of the form 
factor subtraction, such that there is no dependence on C|J left, i.e. the primitive kernel 
nil = 0. 

The "tadpole" contractions shown in Figure 3 (c) and (d) can be written as an 
effective flavour-changing two-point vertex 



V2 



pb 



p=u,c 



{C5 + N,CQ)AmlD{l + j,)b, 



(33) 



where A is a dimensionless, ultraviolet- divergent constant. However, the twelve tree 
diagrams with this vertex insertion vanish in their sum, so there is no contribution from 
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this class of diagrams. We also find that the diagram classes (b) and (f) in Figure 3 
actually vanish. 

Penguin contractions to the hard-spectator scattering kernels relevant to the QCD 
penguin amplitude have previously been calculated in [7]. In this paper the diagram 
classes (b), (c), (d) and (f) in Figure 3 are not considered (fortunately they all vanish), 
but diagrams (a5), (a6), (e3), (e4) and (e5) are also omitted, which is incorrect. The 
correct procedure includes these diagrams together with the form-factor subtraction 
discussed above. This seems to lead to significant numerical differences between our 
results and those of [7]. 

Electroweak penguin contractions (Figure 1 (f) — (m)) 

Electroweak penguin contributions from topologies specified in the first two lines in 
Figure 3 can be obtained by straightforward adjustment of colour factors from the cor- 
responding QCD penguin diagrams. 

The evaluation of the set of diagrams shown in the third line of the figure, which 
contributes to af^EW^ requires some explanation, since the photon has small virtuality. 
When the vector meson M2 is transverse, this leads to an enhancement [11], but for 
the longitudinal polarization state considered here the photon propagator cancels, and 
the diagrams contribute to the matching coefficients of SCETi four-quark operators as 
shown for the non-spectator kernels in [9]. 

We calculate this contribution by keeping the up-going quark and anti-quark on-shell 
{Qi = Q'i ~ 0)' t>ut we assign them a small transverse momentum q2± = —qi±, leading 
to a non-zero virtuality of the intermediate photon. Then a straightforward calculation 
shows that we may perform the substitution 

^ f^?"^ - (1 - e)^) <l{qihMl2) - -2^ f A - ^) X^X (34) 

in these diagrams. The term contracts to zero with the remainder of the diagrams 
(after summing them all) as required by gauge invariance; in the second term the photon 
pole is manifestly cancelled. Nevertheless, one cannot set = at this point, since the 
loop integrations result in In . These infrared logarithms are subtracted by performing 
a SCET matching calculation. With 7^ the SCET matrix elements do not vanish and 
contain ln(g^/z/^). After matching the hard-scattering kernels are finite as — > 0, but 
depend on the ultraviolet subtraction scale z/, which is related to electromagnetic scale 
dependence of the longitudinal neutral vector meson decay constant. This scale depen- 
dence arises because (contrary to a widely held opinion) the electromagnetic current has 
an anomalous dimension in QED, which is precisely related to the penguin contractions 
as noted in [9]. See also [22] for an explanation of this point. 

3.2 Results for the primitive kernels 

The analytical results for the one-loop kernels are as follows. 
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Vertex contractions (Figure 1 (a),(b)) 

The kernels ri and r2 can be found in [4]. The two new vertex kernels are related to 
these by 

r3{u,v) = -ri{u,v) + {4Cf - CA)(^r2{u,v) + (35) 
6 

ri{u,v) = r2{u,v) + -. (36) 
The colour factors are Cp = (N^ - l)/(2iVc) = 4/3 and Ca = N^ = ^. 
QCD penguin contractions (Figure 1 (c) — (e)) 

All kernels can be expressed in terms of the four functions Si, s^, sf, s'^. With 

F, , 141 1 8{u + uv)s 2{l-2v)s^. _ , 

Ssis) = -f-- - = Co{s, u, v) 

9 u uv iw^v uv 

fSus 6-4mi\ ^ ^ /4 8s 2 \ . 
+ [s^ - j ^) + is + 3^ + ^ j 



, 5 1 1 2s vs ^ . , 
4{s) = -T^- + TT— + — + — Co s, u, v) 
18 u 6UV n^f uv 



we have 



/ 1 2{3-v)s\^ , , /4s l + 2w\ , , 



si{s) = -- — ^ _^^^'^ Co(s,M,^) + ^(Go(s,Mi;) - Go(s,m)), (39) 

u 2uv uv uv 

s\{s) = -si(s), (40) 

S2{s) = S;{s)-{cf-^\^, (41) 



2 /3m 



Ca\ 2 



= s,{s) - {ACf - Ca)sM - [Cf - -fj — , (42) 

S3is) = CFS^{s) + CAsi{s), (43) 

4(1)= (2^^-^)4(1), (44) 
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S4(s) = Si{s), 

4(1) = -- - ^^^^^Co(l, u, v) - -(Go(l, uv) - Go(l, u)), 
uv uv uv 



mi = 0. 

Here 



a/4s - X - ie 
Go(s,xj = -= arctan- 



X \jAs — X — ie' 



, , 2 / 2 / ^^1/ 2 / ^ 

Go(s, X, 7/1 = — arctan J- arctan J 

x{l — y) \ \l As — xy — le y 4s — x 

Electroweak penguin contractions (Figure 1 (f) — (m)) 

51, Ew(s) = iVc^sf (s) - 2si(s) - Jr^, 

4,Ew(s) = Sl,Ew(s), 

2 

52, EW(S) = (s) - — , 
^.Ewl'^) = S2,Ew(-5) -4si(s), 

53, EW(-S) = sf (s), 

4,Ew(l) = 24(1), 

54, Ew(-5) = iVc(sf (s) - 2Si(s)), 

■54,Ew(l) — 0' 
"^l.EW = 0, 
"l2,EW = 0, 



S5,EWlSj 



A/'c(6 + 2i7r-21nw-21n--^) (s = 0), 



iVe(6 - 4Go(s, v) + 4sCo(s, t;, 0)) (s ^ 0), 

S8,Ew(s) = S5,Ew(s), 



"58,EWI'5J = -Ss.EWlSj, 
"23,EW = 4A^c, 
'^4,EW = 0. 
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Parameter 


Value /Range 


Parameter 


Value /Range 


. (5) 

A— 

MS 


0.225 






rric 


1.3 ± 0.2 


f^hc 


1.5 ± 0.6 


ms{2 GeV) 


0.09 ±0.02 




0.21 ± 0.02 


{rriu + md)/ms 


0.0826 


U [h] 


0.131 [0.16] 


nil, 


4.8 


/f"(0) 


0.25 ± 0.05 


rhbirhb) 


4.2 


/f^(0) 


0.34 ± 0.05 


\Vcb\ 


0.0415 ±0.0010 




0.32 ±0.05 


\Vub/Vcb\ 


0.09 ± 0.02 


Xb{1 GeV) 


0.35 ± 0.15 


7 


(70 ± 20)° 


ai{l GeV) 


1.5 ± 1 




1.64 ps 


(72(1 GeV) 


3±2 




1.53 ps 


af (2 GeV) 
a^'^(2 GeV) 


0.06 ±0.06 
0.2 ±0.15 



Table 2: List of input parameters. Dimensionful parameters are given in units of 1 GeV. 

The kernels se.EW, sy^Ew, ^yew vanish. The kernels Sgg-^^ and SQ■^^^ vanish after the 
rearrangement of Wilson coefficients discussed after (31). The scale z/ in S5,ew is due to 
the electromagnetic scale dependence of the neutral vector meson decay constant. In the 
following we no longer distinguish u from the matching scale /i in the other expressions. 

4 Numerical penguin amplitudes 

In this section we study the numerical values of the calculated corrections to the penguin 
amplitudes af(MiM2), af^-yvl^i^s) {i = 3,4), and provide updated results for some 
phenomenologically important penguin-to-tree and penguin-to-penguin amplitude ratios. 

4.1 Input parameters 

The calculation of non-leptonic decay amplitudes needs several input parameters, which 
we summarize in Table 2. These are fundamental parameters such as the strong coupling, 
quark masses and CKM parameters; the matching and renormalization scale parame- 
ters fib and /ihc (see following subsection); and hadronic parameters related to decay 
constants, form factors and light-cone distribution amplitudes of the mesons. 

There has been significant progress in the calculation of the first few Gegenbauer 
moments of the pion and kaon light-cone distribution amplitudes from QCD sum rules 
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and lattice QCD [23, 24], which leads to a change in the corresponding entries in the 
table compared to earlier analyses.^ Less is known about the moments of the B meson 
distribution amplitude, 

-^—=[ — </)B+(a;;/i), an{^J^) = Xsif^) [ — /x) In" -. (64) 

Afi^/ij Jo UJ Jo LU UJ 

The logarithmic moments enter only in the 1-loop correction to the hard-coUinear ker- 
nel J in (22). However, is very important, because the hard spectator-scattering 
contribution to the decay amplitudes is directly proportional to 1/Xb- The value of 
Xb adopted here follows [1] and represents a compromise between QCD sum rules and 
models of the B meson distribution amplitude that seem to favour larger values [25] and 
data that favour smaller values [4, 9]. Recent evaluations of the 5 — s> vr form factor at 
= also tend to smaller values [26], and we have therefore adjusted the input to the 
number taken in scenario S4 of [9]. The difference in F^^(O) vs. F^^(O) is consistent 
with the large SU(3) breaking found in [23, 26]. 

Parameters for mesons other than pions and kaons or parameters not given in the 
Table are taken from [9]. 



4.2 Renormalization group improvement and numerical imple- 
mentation 

Our main task is to evaluate the integrals 

1 /"^ 

I = ~ 2^BA/i(Q) _/) dudvH]l'P{u,v;fii,c)^Ahi'^ - v; fii,c)(pM2{u; fihc) (65) 

in (20) with Hmi(1 — v; fi^c) given in (22). Hard spectator scattering depends on the 
two scales fib ~ f^b and /ihc ~ y/rUhA. In the previous equation we indicated the renor- 
malization scale in the arguments of all quantities. To avoid formally large logarithms 
of the ratio of the two scales we need H^l'^ at /^hc in this equation. On the other hand, 
our calculation refers to matching the effective weak Hamiltonian to SCETj at the scale 
Hb, i-e. the scale /i in Section 3.2 is fib- In the following we explain the evolution of 
the hard-scattering coefficient to /ihc, and the subsequent numerical evaluation of the 
integrals /. 

The renormalization-group evolution of the hard-scattering coefficients is the same 
for all coefficients and can be written as 

Hll'P{u,v-fi^,) = e-^^^^-'^'-) [' du'dv'UBL{u,u'-fib,fihc) 

Jo 

xf/||(l -v,l- v';Hb,IJ'hc)Hll'^{u',v';Hb). (66) 



■'Our kaon Gcgcnbauer moments are defined as {K\ . . . |0) matrix elements rather than (0| . . . \ K) and 
therefore the odd moments have opposite sign compared to those papers. 
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This follows because the evolution function factorizes into a term from the Brodsky- 
Lepage kernel [27] for the evolution of light-meson distribution amplitudes, and a Su- 
dakov suppression factor times the evolution kernel for the so-called SCET Bl-type 
currents [5, 6]. (The arguments of U\\ follow the convention of [6], where they refer to 
the gluon momentum fraction in the Bl-type operator, which is 1 — v.) The Ubl factor 
can be used to evolve 0m2(w; fihc) to the scale fib, so (65) reads 

5(M.,A%c) dudv' <Pm, {u; fib)Hl'/{u, v'- fib) 



2F^A./i(0) 

X / dvUi\{l-v,l-v']fib,fihc)'^Mi{l-v]fJ'hc)- (67) 



The hard-scattering coefficients can be divided into two terms, 

H^' = H^tl + H'^'l (68) 

referring to the vertex contractions (primitive kernels r^) and penguin contractions (the 
others). The renormalizat ion-group evolution and numerical evaluation of the first part 
is done exactly as described in Section 4.2 of [4]. In the following we discuss only the 
second part. This part begins at order (or aem), hence it is sufficient to evaluate S^i 
in the tree approximation (order a^) for the hard-collinear kernel J resulting in 



nCp a^(/ihc)/A/i/B( /ihc) 0Mi(t'; /ihc 

Nc mbXsifJ'hc) 



iMi(i -^^;/^hc) = — 77 — ^ = • (69) 



Now we expand (f)Mj^{v; fi^c) into Gegenbauer polynomials 

oo 

<PmAv; /ihc) = 6vv (/ihc)C(3/'n2^ - 1), (70) 

n=0 

where a^^'^{fihc) are the Gegenbauer moments (ao^^(yUhc) = 1), and define 

/■I 

Cn{v] fib, fJ'hc) = dw6wCj^^'^\2w-l)Uii{l-w,l-v]fib,fihc)- (71) 

JO 

Inserting this into (67) we obtain 

m,F^A/,(o)AB(/ihc) 2iV, 

oo I 

xJ^^n^f^hc) dudv(j)M2{u; fib) Hpf^{u,v; fib)Cn{v; fib, fJ'hc)- (72) 

It remains to calculate the Cn{v] fib, f^hc) ■ This can be done by solving numerically the 
integro-differential equation 

fi-^Cn{v;fi,fii,c) = -J^ dwjwil -v,l -w)Cn{w;fi,fiY,c) (73) 
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subject to the initial condition C„(f ; /ihc, A^hc) = GvC^^"^^ {2v — 1), which follows from the 
defining equation (71) and the renormalization-group equation for U\\. Here 7|| (r,r') is 
the leading-order anomalous dimension of the Bl-type currents [5, 6] (Eq. (99) in [6]). 
When this is done, the leading-logarithmic terms are summed to all orders, and (72) is 
free from formally large logarithms of the ratio of and A/m^A. 

We truncate the Gegenbauer expansion of 0Mi and after the second moment, so 
(72) is a sum of nine terms, each proportional to a product of two Gegenbauer moments. 
For fixed /ihc and fib the evolution equation and the remaining integrals can be solved 
numerically. In practice, we wish to construct a code that allows us to vary the scales 
freely in order to estimate theoretical errors. The numerical evaluation is then too time- 
consuming. We therefore proceed as follows. First we solve the evolution equation (73) 
by successive approximation up to the second order. Since the leading-order anomalous 
dimension 7|| depends on fi only through ^^(/i), it is convenient to introduce a variable 
X via 

u— = /5 a, — = (74) 

dfi das 47r dx 

with X = for /i = /ihc. Then we solve (73) up to and including terms of order x"^. Given 
/ife and /ihc, the required value of x is 



X 




as(/ihc) I3q 



(75) 



with /3o,i the first two coefficients of the beta-function, here for four flavours. It would be 
consistent to drop the Pi terms, but since we always use 2-loop running of a, in our code, 
we also include the corresponding terms here. We have checked that the second-order 
approximation to the numerical solution of the evolution equation always provides an 
adequate approximation, even for the largest scale hierarchy allowed by Table 2, when 
X reaches about —0.06. 

We thus obtain Co. 1,2(^5 x) as second-order polynomials in x, which can be integrated 
numerically coefficient by coefficient in (72). A final complication arises due to the 
dependence of some of the Hp'^{u,v] fib) on the charm quark mass, more precisely on 
rric/mb, which should also be allowed to be variable. We deal with this by integrating 
(72) for several values of mc/rrib, and then generate a polynomial fit which approxi- 
mates the dependence on mc/mb in the relevant interval. Thus the final result for / 
is represented as a sum of nine terms, corresponding to the coefficients of products of 
Gegenbauer moments. Each term is a second-order polynomial in x, and in some cases 
also a polynomial in rric/mb with numerical coefficients, such that the dependences on 
all parameters can be rapidly evaluated. 

The numerical results below include the renormalization-group evolution as described 
here. We have, however, found that the effect of summing logarithms is not very im- 
portant. A simpler implementation leaving out evolution would not lead to essential 
differences in the sense that any difference to the evolved results is smaller than other 
uncertainties. 
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4.3 Amplitudes 

The final result for the af parameters in (20), now including all contributions, can 
be written in a notation similar to Eq. (35) of [9], which makes explicit the contribu- 
tions from the form-factor term in the factorization formula (1), from hard spectator- 
scattering, separately from tree- and 1-loop corrections, and from penguin contractions. 
To simplify the notation we use the identification ai = ai, 02 = 02, = afi, 04 = af^, 
= afjij, Qj = af/jEW' ^9 — "^slew; '^10 = '^4Lew- We also give the coefficients 
Og 8 related to the power-suppressed "scalar" penguin amplitudes for completeness, even 
though the 1-loop spectator-scattering correction has not yet been computed, and add 
the tree-level twist-3 spectator scattering term to conform with the conventions of [9]. 
The general formula is 



af (M1M2) = Q 



Ci±i , Ci±i as 



+ 



V,{M2) + P[{M2) 



+ 



iV2 mbF^A^i(0)AB 



+ Ci HVl^'^ (M1M2)) + i/Pf (M1M2) 



(76) 



(Recall that F^*^i(0) = /+*'^H0) for pseudoscalar Mi and A^*^^(0) for Mi = V.) The 
upper signs apply when i is odd, the lower ones when i is even. The non-spectator terms 
Vi(M2), Pf(M2) are listed in [9] (see also [2]). The tree-level spectator-scattering term 
Hi{MiM2) differs from the convention of [9] by an overall factor. Introducing 



Am = 
we have 



3x 



n „M 



a 



n ' 



A 



M 



n=l 



3x 



H,{M,M2 







1 + 

n=l 

1-4,9,10 

5,7 

6,8. 



(77) 



(78) 



Here Xh models a non-factorizable power correction (since ~ l/^rife) as defined in [9] 
[Eqs. (50) and (63)]. The new 1-loop spectator-scattering correction is contained in the 
objects if\^/^^(MiM2), HVP{MiM2), HPi{MiM2), which we now summarize in terms of 
the previously defined Hjl [Eq. (24-31)] and primitive kernels. The vertex contractions 
(primitive kernels r^) contribute 



HV^^\MiM2) 



Pi(MiM2) + Am,J(Mi) 

P3(MiM2)-Am,J(Mi) 

unknown 



1-4,9,10 

5,7 

6,8, 
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HVj:'^\MiM2) 



R2iMiM2) 

i?4(MiM2) 

unknown 



1-4,9,10 



i = 5,7 
i = 6,8. 



(80) 



Here Ri{M\M2) denotes the integrated vertex-contraction kernels 

1 T in v] 

Ri{MiM2) = - / dudv(l)MAv)(l)M2iu) - ■ 



9 Jo 



(81) 



and the next-to- leading order hard-coUinear correction, J (Mi), is given in Eq. (50) of [4]. 
The above expression for Ri{MiM2) does not include the renormalization-group sum- 
mation of logarithms, since the final formula is complicated. The contributions from 
the penguin contractions to spectator scattering are simply (72) up to a normalization 
factor, 



HP[{M,M2) 



^ n=0 



1 

X / dudv (f)M2 (u; fib) 



Nc rrU,p 



Hpfk{u,v;fib) 



C„(w;/ife,/ihc), 



(82) 



which includes the log-summation. Again, this result does not apply to the power- 
suppressed penguin amplitudes, i = 6,8, for which the radiative corrections are currently 
not known. 



4.3.1 The tree amplitudes 0:1,2 



The "tree" amplitudes ai,2(7rvr) have already been discussed in [4]. With our up-dated 
input parameters, their values read 



ai(7r7r) = 1.015 + [0.025 + 0.012i] 



V 



sp 



0.485 



a2[7nT) 



{[0.020]lo + [0.034 + 0.029i]Hv + [0.012]tw3} 
0.975t°:°?^ + (-0.017tnf ^ 
0.184- [0.153 + 0.077i]y 



i.05lJ^5 



(83) 



+ 



' sp 



0.485 J 



{[0.122]lo + [0.050 + 0.053t]Hv + [0.071]tw3} 



0.275 



+0.228 
-0.135 



(-0.024;°; 



081 ) 



(84) 



In these expressions we separated the tree {a^, first number), vertex correction {as, 
indexed by V) and the spectator-scattering correction (remainder). The latter is further 
divided into the tree (a^, indexed LO), 1-loop (a^, indexed HV), and twist-3 power 
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correction (the Xh term in (78)). The theoretical error in the last line of each expression 
is computed from the ranges in Table 2 added in quadrature together with the error from 
Xh, which is treated as described in [9]. The most important parameter uncertainty is 
encoded in the combination 

which normalizes the spectator-scattering term as can be seen from (76). 

The dependence on the final state mesons is rather small. For instance, the difference 
between ai^2{'^K) and ai 2(7r7r), often called 'non-f actor izable' SU(3) breaking, is about 
10% for 02 and the imaginary part of Oi (and much smaller for the real part of Oi), 
since the recent estimates of the first Gegenbauer moment of the kaon light-cone distri- 
bution amplitude favour small values. The dominant SU(3)-breaking effect can thus be 
estimated from the above expressions by the dependence of r^^ on fu^lF^'^^^ifi). 

4.3.2 The QCD penguin amplitude cx^ 

The QCD penguin amplitudes a^"'^ receive a 1-loop spectator-scattering correction due 
to the primitive kernels ri, r2 and the newly computed Si . . . S4. We find 

al{'KTi) = -0.029 - [0.002 + 0.001z]y + [0.003 - 0.013z]p 



+ 



485 J {[°-°°^]lo + [0.001 + O.OOU]hv - [0.000 + O.OOU]hp + [O.OOlJtws} 
Omit^o'ot + {-O.Ol2t'oZ)^ (86) 



aKmr) = -0.029 - [0.002 + 0.001z]y - [0.001 + 0.007z]p 



+ 



485J {[°-°°^]lo + [0.001 + OmUUv + [0.000 - O.OOO^Ihp + [0.001]tw3} 

= -omst'oZl + (-0.006^°:°°^^ (87) 

Here "P" and ^^HP" denote the contributions from penguin-contraction diagrams to the 
form-factor and spectator-scattering term in the factorization formula. Numerically the 
new corrections, labeled ^^HV" and ^^HP" , respectively, are very small. ^ As the large 
Wilson coefficient Ci is involved in the new penguin correction (the term CiSs{sp) in 
(24)), this result is somewhat surprising. Closer inspection shows that 

[<(vrvr)]HP = [0.0074 + 0.0060i]c^Ci - [0.0073 + 0.0053i]cpCi - 0.0002 + O.OOOU, 

[al{n7r)]HP = [0.0039 - 0.0000^]^^^ - [0.0034 - 0.0006z]c^Ci - 0.0004 - 0.00022, 

where the terms labeled "C^Ci" and "CpCi" arise from the term CiS^lsp) separated 
according to the two colour structures Ca and Cp (primitive kernels s^{sp) and (sp)). 



^Thc terms labelled 'P' and contain the charm penguin contractions. There is no evidence 

from this calculation that these effects are large. 
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Figure 4: LO, NLO and NLOgp value of 04 in the complex plane, has a sizeable phase. 
The NLOsp point includes a theoretical error estimate. 



Thus, there is an almost complete cancellation between the two terms, which individually 
would have resulted in a large correction to the QCD penguin amplitude. Inspecting (37), 
(38) some terms appear with the small colour factor Cp — C*a/2, but it is unclear to us 
whether there is an explanation for the near-completeness of the numerical cancellation 
or whether it is accidental. 

A graphical representation of 1-loop corrections to 04 is given in Figure 4. The 
leading-order (LO) point corresponds to naive factorization. The point labelled "NLO" 
includes all as corrections. The point including the error bars is the result of our cal- 
culation that adds the 1-loop spectator-scattering correction. We refer to this partial 
next-to-next-to-leading order (NNLO) calculation as "NLOgp" , since the NNLO term is 
a 1-loop correction in the spectator-scattering sector. The point "G" corresponds to a 
choice of inputs made in [4] to achieve an improved agreement with the B ^ tttt decay 
rates. This input parameter set uses = 200 MeV and 03 = 0.3 near the boundaries 
of the assumed regions instead of the central values given in Table 2, which leads to an 
increase of the spectator-scattering contribution. 

The figure shows that the new correction has a small effect on the real part, but can 
be relevant to the imaginary part depending on parameter values. As was found in [4] 
for the parameters ai, 02, perturbation theory is well-behaved. These observations apply 
to all final states, in particular, "non-factorizable" SU(3) breaking turns out to be quite 
small, as discussed above. 

To gauge the impact of the spectator-scattering correction on branching fractions 
and CP asymmetries, we recall that 04 always appears in conjunction with the power- 
suppressed penguin amphtude r^Og. Since we have not computed 1-loop corrections to 
the short-distance coefficients of power-suppressed operators, the value of r^ag differs 
from previous analyses only due to the updated input parameters. Furthermore, since 
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the tree-level spectator scattering contribution vanishes for ag, the following numbers are 
solely due to the naive-factorization and 1-loop contributions to the form-factor term in 
the factorization formula. We have 



TTTT 



TTTT 



-0.055 - [0.001]y - [0.004 + 0.020i]p 

-0.060^°:°}^ + (-0.0201°:°°^)^, 



(88) 



-0.055 - [0.001]y - [0.009 + 0.010z]p 
-0.065+°-°i2 



l+0.004\ 



_0.oi9 -r (-0.010™)2. (89) 

It is worth noting that r^ag(7r7r) is numerically larger than 04(717?). On the other 
hand, rJag(MiM2) is strongly suppressed when M2 is a vector meson, where the naive- 
factorization contribution is zero. We shall discuss some of the phenomenology associated 
with the penguin amplitudes in Section 5. 



4.3.3 The flavour-singlet QCD penguin amplitude 

The flavour- singlet QCD penguin amplitude is not relevant to the 7r7r final state. We 
therefore give numerical values for 



a^iKcj)) = 0.001 + [0.005 + 0.002z] 



V 



' sp 



0.435 



{[0.003]lo + [0.001 + 0.001i]Hv + [0.002]tw3} 



001 

U.UUi„o.005 



l-0.002^ 



af(J?0) 



(O.OOllo.oos/^ 
0.005 - [0.002 + 0.003i]v 



(90) 



+ 



sp 



0.435 



{[0.004]lo + [0.002 + 0m2i]Hv + [0.002]t„3} 



001+°-°°'^ + 

U.UUi_Q QQ4 -t- 



-0.000 



h0.005\ 
-0.003 J 



(91) 



valid for p = u and p = c. Here the 1-loop spectator scattering term ^^HV" is numerically 
of the same size as the other terms, and gives a large effect. This is typical for corrections 
to the smaller colour-suppressed amplitudes. We recall from Section 2 that the calcu- 
lation of the flavour- singlet amplitudes afg is not complete when M2 is a pseudoscalar 
meson, since we did not compute the matching coefficients of two-gluon operators. Even 
more important conceptually, there is another term in the QCD factorization formula 
which does not factorize in the form of (1), which has been estimated to be of similar 
size as the individual terms in (90), (91) [12], but which is very uncertain. 

Because the flavour- singlet amplitude = as =F 05 is so small, it does not play 
an important role in the phenomenology of the more prominent charmless final states, 
in particular in the explanation of the large rj'K branching fractions [12]. The large 
spectator-scattering correction found above affects final states such as r]'K* or such with 
very small branching fractions, where the leading contribution from is suppressed. 
However, the predictions for such decays carry large theoretical uncertainties. 
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4.3.4 Electroweak penguin amplitudes 

The electroweak (EW) penguin amplitudes are small, yet they contribute significantly to 
isospin-breaking effects in decays to ttK final states. The "colour- allowed" EW penguin 
amplitude is af ew = ^^9 -F c^? (upper sign for tiK). We find 

a^{K?c) = 10^2 X (0.009 + [0.002 + 0.003i]v + [0.038]p 



sp 



0.435 



{[0.005]lo + [0.003 + 0.004i]w - [0.020 + 0.016i]^^p + [0.003]tw3} 



10-^ X (0.0581°:°?^ + (0.0151^ 



-+0.010N 



0Q6J'' h 



(92) 



a^jiKn) = 10"^ X 0.009 + [0.002 + 0.003i]y + [O.Olljp 



sp 



0.435 



{[0.005]lo + [0.003 + Om^HV + [0.002 - Omii^P + [0.003]tw3} 



10"^ X (o.oioi[J:[]}^ + (o.ooollj 



+0.003 >| 



006 J '' h 



(93) 



a^iKn) = 10~^ X - 0.909 - [0.023 + 0.011i]y + [0.038]p 



+ 



sp 



{[0.017]lo + [0.028 + 0.025i]Hv + [0.020 + O.OWijHP + [0.010]tw3} 



0.435 

10-^ X ( - 0.819l°:°«° + (0.029l°:°|)z 



(94) 



a'giKTr) = 10'^ X - 0.909 - [0.023 + O.OIH]^ + [0.011]p 



+ 



sp 



{[0.017]lo + [0.028 + 0.025i]Hv - [0.002 - 0.001i]HP + [0.010]tw3} 



0.435 

10"^ X ( - 0.868t[!:[J| + (0.015t 



+0.043\ 

o.oisJ 



(95) 



The penguin-contraction term "i/P" originates from the diagrams shown in the third 
line of Figure 3. Both "P" and "iJP" coincide for and ag, such that they exactly 
cancel in the combination ag — ay appropriate to final states with pseudoscalar M2, but 
add up for vector M2. The 1-loop spectator scattering correction is quite important 
for Oy, and a^ is seen to even be dominated by the penguin contractions. However, 
what is relevant to decay amplitudes is a^^EW: and therein the 1-loop corrections are 
overwhelmed by the large tree contribution to ag. 
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Figure 5: The electroweak penguin amplitude Oiq receives a large correction. The am- 
plitude is almost identical for p = u and p = c, therefore only the case p = u is shown. 

The colour-suppressed EW penguin amplitude a^-^w = '^lo i^x'^s consists of afo and 
a power-suppressed amplitude r^^ag, which is the electroweak equivalent to r^ag. We find 

r^a^inK) = 10'^ x (0.063 + [0.000]v. + [0.004 - 0.016z]p ) 



10- X (0.0681°;°^^ + (-0.016™)z 



+0.010\ 



(96) 



r^aliirK) = 10"^ x 0.063 + [0.000]v^ + [0.002 - 0.010i]p 



10-^ X (OM5t'ofi 



(97) 



a'^Q^nK) = 10~2 X - 0.161 + [0.135 + 0.073z]y + [0.026 - 0.01H]p 
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alo{nK) = 10^2 x - 0.161 + [0.135 + 0.073z]y + [0.023 - 0.006«]p 
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The coefficient is similar to the colour-suppressed tree amplitude a2- The non- 
spectator term (ffist line in each expression) is strongly suppressed due to a cancellation 
between the naive-factorization term and the 1-loop correction, which is large due to 
the absence of colour-suppression. The final number is largely from spectator-scattering 
which obtains significant 1-loop corrections. The result is displayed graphically in Fig- 
ure 5. 

5 Amplitude ratios 

With the improved parameters we are in a position to calculate complex ratios of 
strong amplitudes such as P/T or C/T. In this section we discuss a few prominent 
examples. 

5.1 The PP, PV, VP, and [VV]l QCD penguin amplitude 

We recall that in physical decay amplitudes the parameters a^, af, and the penguin 
annihilation amplitude /^f always appear in the same linear combination [9] 

al{MiM2) = a4(MiM2) ± r^'ae{M,M2) + /3f (M1M2), (100) 

where the upper (lower) sign corresponds to the case Mi = P {V), and cannot individu- 
ally be confronted with experiment. The annihilation contribution (3^ is incalculable in 
factorization, since it contains an endpoint divergence; we therefore resort to the model 
defined in [2], parameterizing it in terms of a complex parameter 

X^ = ln^(l + f?^e^*-) (101) 

with Ah = 500 MeV a hadronic scale and f?A = by default. An additional theoretical 
error is assigned to any observable by setting qa = ^ and allowing the phase 0a to take 
arbitrary values. It turns out that this is almost always by far the largest theoretical 
uncertainty for direct CP asymmetries, and among the largest uncertainties for branching 
fractions of penguin-dominated decays, especially those with vector mesons in the final 
state. The expressions for the annihilation amplitudes are given in [9] and, for VV final 
states, in [28]. 

To assess the validity of the QCD factorization framework, we show in Figure 6 
the ratios a^(MiM2)/_(ai(7r7r) + a2(vr7r)) (M1M2 = ttK, ttK*) and al{MiM2)/{ai{pp) + 
a2{pp)) (M1M2 = pK, pK*), the ffist two of which have been previously considered 
in [9]. (For pK* and pp, only the longitudinal polarization amplitude is considered 
in the following.) The result of the calculation is represented by the dark point with 
error bars. The nearly circular contours around this point show the variation of the 
theoretical prediction when the phase of the annihilation model is varied from to 
271 for fixed qa = 1,2,3 (inner to outer circles). The blue square corresponds to the 
parameter set G, which is defined by = 200 MeV, = 0.3, = 80 MeV, Qa = ^ 
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Figure 6: Comparing the PP, PV, VP, and VV penguin amplitudes to data. The figures 
in the upper row show theory predictions for al{MiM2) / {ai{mi) + a2(vr7r)); those in the 
lower row show dl{MiM2) / {ai{pp) + a2{pp)). Where available, the ranges for modulus 
and phase extracted from data are also depicted. See text for explanations. 
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and = -55° (PP), 0a = -20° (PV), (pA = -70° {VP), qa = {VV), following the 
favoured parameter set S4 of [9]. 

As far as data is available, the amplitude ratios can be obtained with little theory 
input from the well measured (CP-averaged, longitudinal) branching fractions Br(P~ — > 
7r-i?(*)0), h{p-K*^) ■ Br(P- ^ p-R*^), Bt{B~ ^ rr-n^), hip-p') ■ Br(P- ^ p~p^), 
and the rate and direct CP asymmetry in ■k~^K^*'^~ (cf. [9], (75) and (77)). The 

darker rings are due to the experimental errors in the branching fractions and the lighter 
ones include also the uncertainty of \Vub\- The angles of the wedges involve the CP 
asymmetry measurements (darker region) while the lighter region also includes the error 
on 7/ The wedges opening to the right are ruled out (or at least disfavoured) by the 
fact that the measured values of the Fleischer- Mannel ratios [29] 

R^K_ nB'^n^K-) _ r(P°^vr-^JC") _ 

(102) 

are both less than unity. 

Three messages can be read off from the figure. (1) The magnitudes and phases pre- 
dicted for the amplitude ratios agree reasonably well with data, as indicated by the error 
bars and the small onion-shaped regions. A large annihilation amplitude is disfavoured, 
since it would require fine-tuning of the phase to satisfy the experimental constraints, 
but some annihilation contribution appears to be required, especially for the PV ampli- 
tude. (The apparent smallness of the annihilation uncertainty in the pK* plot is due to 
cancellations in the crude annihilation model employed here.) There is a tendency of the 
predicted magnitude of the penguin amplitude without annihilation to be smaller than 
the data. This difference is about 0.02 to 0.03 independent of the spins of the final-state 
mesons, except for the VV final state where no such difference exists. (2) The magni- 
tude of the penguin amplitude is predicted much smaller in the case containing a vector 
meson in the final state, either due to the smallness of Og {nK*, pK*) or a cancellation 
of 04 and ag (pK). This is also reflected by the data. (3) A priori the ratio could have 
lain anywhere in the complex plane shown in the figure. The agreement found hence 
constitutes a highly non-trivial check of the qualitative and quantitative predictions for 
penguin amplitudes in the factorization framework. In this graphical representation the 
well-known difficulty to account for the small direct CP asymmetry in 7i~^K~ in 

QCD factorization is seen as a small offset of the theoretical calculation from the left 
wedge in the first panel of the figure. 



5.2 P/T and C/T 

Penguin-to-tree and other ratios are also of phenomenological importance in P — irn, 
B —>■ Tip, B — > pp, and other decays. The ratio P^tt/Ptttt, for instance, can be determined 
solely from the time-dependent CP asymmetry in P ^ vr+vr" for given values of the well 

^In calculating the wedges, an additional, small theory uncertainty on the ratio Q;i(7rA''^*)'')/(ai(7r7r) + 
a2(7r7r)) is not included. 
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Ratio Value/Range 



Value G 
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n 1 Q4-^-0.199 



;o.o29i°:l^^)z 

-0.024t°:l??)z 



0.691 + 0.165i 
0.344 + 0.042Z 
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T, 



pTV 



T 



0.869tJ]: 



+0.275 
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0.014lg; 



058 N 
057 J 



0.945 - 0.004i 
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Table 3: Amplitude ratios for the nir, pp and np final states. In the case of pp the 
ratios of longitudinal polarization amplitudes are given. The third column gives the 
preferred-parameter-set G value. 

measured mixing phase (pd, equal to 2/5 in the Standard Model, and the CKM angle 7. 
Conversely theoretical predictions for these ratios allow to extract the angle 7. Similar 
relations hold for final states ir^p"^, p^p~, 7i~^K^*^~, etc. Here we present numerical 
values for a number of these ratios. 

In such phenomenological studies it is convenient to define the colour-allowed tree 
amplitude T, the colour-suppressed tree amplitude C, and the penguin amplitude P as 
the hadronic amplitudes multiplying the different CKM structures in the decay ampli- 
tude. In this convention, the name of an amplitude derives from its leading contribution, 
such that T ~ ai, C ~ 02, and P ~ a^, but sub-leading terms can make an important 
difference. For the following discussion of -B — > nn, pp, np decays, we define T, C, P 
through 

A-^o^^^^_ oc V*aVubT^p + V^aVcbPnp, 
-2^B"^^OpO oc V*dVub[C^p + Cpn] + V*aVcb[. ■ •]• 



(103) 
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These definitions are related to the af (and annihilation) amphtudes by comparing the 
above equations to [9], (A. 13). The proportionahty factor is the same in all three lines 
and therefore irrelevant, since we shall only consider ratios. Analogous definitions apply 
to utt and pp. 

Our results for the amplitude ratios are given in Table 3. The first column of num- 
bers provides the default result with errors, the second column the value in parameter 
set G that we currently consider as our "best" result. The theoretical calculation has 
significant uncertainties. The dominant error of P/T always arises from the weak annihi- 
lation model and cannot be reduced by further calculations. The error due to neglected 
higher-order corrections is difficult to estimate, but there may be a sizable shift of the 
central values of P/T due to the uncalculated O^a"^) correction to the power-suppressed 
"scalar penguin" amplitude Cg. That correction, like discussed above, involves the 
large Wilson coefficient Ci at one loop, but there may not be a numerical cancellation 
analogous to the one observed for 04 in Section 4.3.2. For C/T, insufficiently known 
input parameters such as and , and to a lesser extent the twist-3 correction Xh, 
are responsible for the bulk of the uncertainty. 

There is a clear hierarchy of the penguin-to-tree ratios that has already been discussed 
in [9] and in the previous subsection. The difference in C/T for the various final states is 
in fact a reflection of the same hierarchy of penguin amplitudes. While 01^2 are roughly 
the same for vrvr, vrp, pp (longitudinal polarization), the amplitudes T,C contain the 
up-penguin amplitude, such that T = ai + + . . ., C = 02 — + . . .. It is worth 
noting that C/T can acquire a signiflcant imaginary part from the 1-loop vertex and 
spectator-scattering correction, and the large phase of the up-penguin amplitude. It will 
be interesting to detect these characteristic features in the experimental data. 



5.3 Relating electroweak penguin to tree amplitudes 

The electroweak penguin amplitudes afEw '^4EW related to the tree amplitudes 
Q!i_2 under certain assumptions [30, 31, 32]. Therefore, more than the calculated values 
of the electroweak amplitudes, the deviations from these relations are of interest. The 
assumptions made in deriving the relations are the neglect of the electroweak penguin 
operators Qr.s, since they have small Wilson coefficient compared to Qg.io at the scales of 
interest, and SU(3) fiavour symmetry, when a relation between tttt and ttK final states 
is involved. In addition there are further assumptions, which amount to neglecting 
the charm and bottom content of the operators Qg^io, and to not considering penguin 
contractions. Since our calculation does not make use of any of these assumptions, it 
is interesting to investigate how well these widely used amplitude relations are satisfied. 
We confine ourselves to the discussion of the -kK and vrvr final states. 
The most solid relation is [30] 



5i 



EW 



ai(vrvr) + a2(7rvr) 2 



cb 



ub 



C9 + C 



10 



C1 + C2 



(104) 



with R^K = fnF^ (0)/(/i^F^^(0)) and a similar relation with vrvr in the numerator and 
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R-kK 1- Let T^NR be the ratio of the middle expression to the right-hand side. We find 

7^NR = (1.02^0:^^)6^(1^1)° (105) 

The closeness to 1 involves a cancellation between an SU(3) breaking on the order of 10% 
and an error of similar size due to neglecting the electroweak penguin operators Qj^s (for 
TTTT we would find 7^NR ~ (0.92tg;[]^) e^(i±i)°). The effect of penguin contractions is below 
the (1 — 2)% level. Note also that the ratio remains real to an excellent approximation. 
The relation holds to much better accuracy than at the time when [2] was written, where 
it was studied without the 1-loop spectator scattering correction. The difference with 
respect to [2] is almost exclusively due to the change in the hadronic parameter ratio 
RnK, which is also responsible for the bulk of the error. 

A second relation follows from SU(3) symmetry [31], which also involves annihilation 
amplitudes. Neglecting such power-suppressed amplitudes, we obtain 



(106) 



ai{Tm) — a2{Trn) C1 — C2 

Let 7?.GPY be the ratio of the left-hand to the right-hand side. We find 

7^GPY = (1.29li)e^(-l-)° (107) 

The deviation from 1, which is bigger than in T^nr, comes about because now the SU(3)- 
breaking correction (for tttt we would find TZq^y ~ 1.13) and the effect of neglecting Ci^ 
add up. Individually, both effects are nearly twice as large as in the first relation, which 
is also refiected in the larger error. Nevertheless, T^gpy remains real to an excellent 
approximation. 

The previous two relations (104), (106) hold under the stated assumptions irrespective 
of the values of the Wilson coefficients Cg^io. Observing that 

- - -0.00896 = (108) 

at /ifo = 4.8 GeV and neglecting the 1% difference between the first three items in this 
expression, stronger relations are obtained in [31], 

0^3,EW = ^^l-: '^4,EW = '^O^Z, (109) 

where due to the SU(3)-symmetry assumption the arguments M1M2 of the Ui can be 
any out of tttt, ttK, Kir and KK, KK. The first of these relations is extraordinarily 
well respected, 

"^'^^^""^^ =(1.00 ±0.01) e-(°-°)°, (110) 

because it involves the colour-allowed amplitudes, and the effect of Qt^^ is very small. 
This holds for any of the above final states vrvr, ttK, KK. In contrast, the second relation 
between the colour- suppressed amplitudes is poor. We find 

"^'^^^^^^ =(0.64+°-^^)e-(^-")°. (Ill) 
fi;a2(7r7r) 
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Our theoretical expectations should be useful to estimate the error incurred by employing 
these relations in data-driven approaches to non-leptonic decays. 

6 Conclusion 

This paper, together with [4], completes the calculation of 1-loop spectator-scattering 
corrections to all flavour-nonsinglet, leading-power decay amplitudes in non-leptonic B 
decays. The calculation shows that factorization works technically: the infrared sin- 
gularities cancel in the matching calculation, and the convolution integrals converge at 
their endpoints. 

The 1-loop corrections are numerically well-behaved. As a general rule, we find sig- 
nificant effects for the colour-suppressed amplitudes 0:2, af and a4E-W' small effects 
for the others. In case of the QCD penguin amplitude a\ this conclusion is reached 
only as the result of a numerical cancellation in the terms proportional to the large 
Wilson coefficient Ci. As a consequence, there is little impact of the newly calculated 
spectator-scattering corrections to penguin amplitudes on the phenomenology of non- 
leptonic branching fractions and CP asymmetries. Further improvement of the calcu- 
lation now requires the calculation of the 2-loop vertex corrections to the form-factor 
term in the factorization formula [33], and an understanding of the power-suppressed 
but large penguin amplitude Og. 

In the final section of the paper we investigated a few amplitude ratios that play an 
important role in determining 7 (a) from tree-dominated h ^ d transitions. Vice versa, 
assuming a value of 7, these ratios may be determined from data and compared to the 
theoretical calculation. A detailed discussion of branching fractions and CP asymmetries 
with account of 1-loop spectator scattering will be presented elsewhere. 
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